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J> ■ Abstract 

o 

^ . 

■ Polarized ferrofluids, lipid monolayers and magnetic bubbles form 

00 
O 

' domains with deformable boundaries. Stability analysis of these do- 



mains depends on a family of nontrivial integrals. We present a closed 

form evaluation of these integrals as a combination of Legendre func- 

Q , tions. This result allows exact and explicit formulae for stability 

O ' 

thresholds and growth rates of individual modes. We also evaluate 

• i-^ - 

X 

' asymptotic behavior in several interesting limits. 
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Stability analysis predicts the evolution of cylindrical polarized domains 
in many systems, including magnetic bubbles ferrofluid labyrinths [0, Q 
and amphiphilic Langmuir monolayers In all these contexts the same 

family of nontrivial integrals arises. We evaluate these integrals in closed 
form in terms of Legendre functions. Our closed form solution allows explicit 
analytic expressions for stability thresholds and individual mode growth and 
decay rates. We also present asymptotic behavior of the integrals in inter- 
esting limits such as large aspect ratio, or high spatial frequency. 

Let us briefly describe one particular physical system. Consider an ini- 
tially circular (radius R) droplet of ferrofluid trapped between two parallel 
glass plates, separated by a distance h. A convenient measure of the droplet 
shape is given by the aspect ratio 

P=^. (1) 

Apply a uniform magnetic field perpendicular to the plates. Varying the 
magnetic field alters the balance between magnetic interactions and surface 
forces. We measure the relative strength of magnetic and surface forces 



1 



through the dimensionless magnetic Bond number 



where M is the magnetization and 7 the hne tension (essentially h times the 
surface tension). 

Jackson, Goldstein and Cebers [Q] calculated the growth (or decay) rate 
of small time-dependent perturbations C, in the radius R of the initially cir- 
cular droplet. They model the global fluid flow and energy dissipation using 
Darcy's law. Let ({(p,t) ~ exp{ant)cos{n(f)), where represents the polar 
angle (figure 1). The growth rate (T„ of the perturbation by mode n is 0] 



Vn{p), which depends on mode number n and aspect ratio p but not on 
the magnetization, is the complicated integral we will shortly analyze. 1] is 
the ferrofluid viscosity. An alternate djTiamical model 0, in which energy 
dissipates locally on the boundary, involves the same integral T>n{p). 
Now we grapple with the integral. Write 




(3) 
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where 

Ux) = d.. X = i. (6) 

■''0 Vx + sin^ 

This integral is proportional to the magnetic energy when the circular droplet 
is perturbed by mode n. It is derived from a double integral around the per- 
turbed boundary (ci; = — 0' is the angular separation of points on the 
boundary). In the late sixties A. A. Thiele derived a number of useful 
results on the asymptotic behavior of the integral In{x) in his study of mag- 
netic bubble stability. Nevertheless, only a systematic expansion of /n(x) as 
a function of the small parameter x has been obtained, an indirect approach 
that requires the use of recursion relations and power series expansions. Even 
more recent publications [2-5] use a similar approach. It would be useful to 
find an explicit closed form solution for In{x). 

To evaluate the integral /„(x) in equation (5), substitute sin^(nu;) = — 
cos(2nci;)] and change variables to 6 = 2uj. The utility of the substitution is 
evident from inspecting the integral representation for the Legendre function 
of the second kind ^ 

1 cos TiO 1 

Qn-iiz) = -= / , de, 3?(n) >--,z> 1. (6) 
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We find a simple and compact expression for tlie solution to the integral 



L(x) 



1 r 



g_i(2a; + l)-Q„_i(2a; + i; 



(7) 



To get closed form expression for Vn{p) we need evaluate the integrals 
In{x) at X = 0. To accomplish this, we expand (5„_i (2x + 1) near x = 



Q„ .(2i+l) 



X 



X+l 



+ 



;i-c)--io^ 



X 



x + l 



il) [n + 



X 



0(x' 



(8) 



where Euler's psi function ip is the logarithmic derivative of the Gamma 
function and C is Euler's constant. Only the lowest order term contributes 
to /n(0), but we will need the first order term in our later discussion. Thus 



4(0) 



(9) 



and finally we rewrite equation 



2 ' J)2 2 p2 



.. (10) 



This is our central result, yielding exact explicit dependence of mode growth 



rates on Bond number Nb-, mode number n and aspect ratio p. 



Figure 2 illustrates the mode growth rate as a function of magnetic Bond 
number Nb for several mode numbers n. We plot the dimensionless growth 
rate 



with (T„ defined as in equation (3). The aspect ratio is held fixed at p = 20. 
Note that mode n = is meaningless because of the constraint of constant 
volume, while mode n = 1 is just a spatial translation, so we consider only 
n > 2. Mode n goes unstable when its growth rate cr„ vanishes, at critical 
Bond number 



The first mode to go unstable as Bond number increases is n = 2. However, 
if the magnetic Bond number is suddenly increased to a value beyond the 
n — 2 threshold, the fastest growing mode n* may have a value greater than 
2. 

In figure 3 we plot the behavior of in terms of the ratio n/p. This 
ratio is proportional to the spatial wave vector k of the perturbation. In 




(11) 
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particular, the wave vector k = n/R and we find 

^ = ^. (13) 

It appears in figure 3 that all curves {2 < n < 15) converge to a limiting 
function as the value of n increases. For large n and p, becomes a func- 
tion only of the ratio n/p. Below we analyze in more detail the asymptotic 
behavior for large n and p. 

In many practical situations, the aspect ratio p or the mode number n 
assume large values [2-5], so we calculate asymptotic expansions taking ad- 
vantage of the closed form expression for T'n(p) derived above. First consider 
large n and p, for fixed ratio n/p. Use the asymptotic limit ^ 

ilmQ„_.(cosh(^)) = iro(g), (14) 
where Koi^q) is the modified Bessel function of order zero, and set 

g = ncosh-^('l + ^U^ + of^V (15) 



Note that 



to get 



lim ip^n + -) ^ log(n), (16) 



T>M^il<oJ'^+Kj?^)-,c]. (17) 



2 1 \p J \ p 



This result is quite general, being valid for large values of p and n, provided 
n/p is kept finite. A related expression may be found in reference 0. 

From equation (12) for the critical Bond number A^^ and asymptotic 
expansion ([17|) for Vnij)) we find, 

Nb ~ / X ^ / X 7- (18) 

Hi«g(f)+^o(?)+c} 

As expected (see figure 3) A"^ is a function only of the ratio n/p. The Bessel 
function Kq simplifies for small or large arguments. In particular, the limit 
of small n/p yields 

i{l-C-log(f)}' ^ ^ 

while in the opposite limit of large n/p. 

Next consider the situation in which p is large, but no assumption is made 
regarding the value of n. This limit is relevant, for example, to the critical 
Bond number of mode = 2, in the limit of large aspect ratio. Typical values 
of p range from the optimal value for magnetic bubbles |^ p = 2 through 



a couple of orders of magnitude for ferrofluid thin films [T^ and reach as 
high as p = 10^ for lipid monolayers ||^. The large p limit can be obtained 
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from equation (^), where the Legendre function of the second kind has been 
expanded to first order in terms of a; = resulting in 



1 f 1 fAn"^ - V 

T^n{p) ~ ^{n^ - 1) log(4p) - 1 - - [^^^ 



(21) 

vahd for arbitrary n. The difference of functions in equation (21) replaces 
the infinite series equation (2.12) in reference [4]. Taking the large n limit 



of (pTD agrees with the small n/p limit of equation ([T7|). 



We conclude by discussing the flat edge limit. Tsebers and Maiorov [|TT 
studied the threshold for instability when a thin layer of a magnetic liquid 
is placed in a homogeneous magnetic field perpendicular to the layer. The 
interface is flat with a small perturbation of wave vector k. To get this flat 
edge limit from our cylindrical geometry, we take radius i? — ^> oo at fixed h. 
Fixed wave vector k implies n and p — oo with fixed n/p as in our asymptotic 
expansion equations (17) and (18). Therefore, replacing n/p with the wave 
vector k as in equation (13) we obtain 

Nb ~ —r --^^ ^. (22) 

2{\og{f)+Koikh) + C} 



This result agrees with the prediction of Tsebers and Maiorov W^ in zero 



gravity. 

Still considering the flat edge limit cited above, we examine the expression 
for the mode growth rate, 



log + Koikh) + C 



- {khy } . (23) 



Taking the limit in which the wavelengths are large compared to the slab 
thickness {kh ^ 1), we expand the Bessel function Kq in equation (23) to 
the second order in kh. In this limit, the critical wave number kc (defined by 
setting a{k) = 0) is 

2 2 

kc ^ ^e^^^e , (24) 
and the fastest growing mode k* (defined by setting da{k)/dk = 0) is 

2 2 2 

ei"^e"^. (25) 

h 

Comparing equations (24) and (25) we see that the ratio of spatial frequencies 
between the critical and the fastest growing mode is 

| = ei (26) 



This result agrees with McConnell [12], who used a simple hydrodynamic 



model to study dynamical stability analysis of a straight edge on a large 
lipid monolayer domain. 
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Figure Captions 



Figure 1: Schematic representation of an initially circular ferrofluid droplet 
of radius R (solid curve) and its perturbation by ( for the n — 5 mode 
(dashed curve). The angular location of points on the boundary is given by 
(j). The droplet volume is kept constant, since the fluid is incompressible. 

Figure 2: Dimensionless growth rate an (see equation (11) in the text), for 
mode numbers 2 < n < 6, as a function of the magnetic Bond number Nb- 
The aspect ratio p — 20. Notice that for each n there exists a critical Bond 
number (when — 0) and for each value of Nb there is a fastest growing 
mode n*. 

Figure 3: Critical Bond number for 2 < n < 15, as a function of 
the ratio n/p. All curves are obtained exactly, combining the closed form 
expression in equation (10) with equation (12). 
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FIGURE 1 
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